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ABSTRACT 



This rsport deals wirh design optimization of Internal- 
Expanding Rim Brakes. A computer program was developed to 
calculate the actuating force, torque, stopping time and 
drum temperature. The drum temperature is calculated by the 
finite difference method. 

A comparison of results has been made using a simplified 
equation that is in common use in engineering texts. 

Numerical optimization is shown to be a convenient tool 
for brake design. 
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SiaSDLS AND ABBREVIATIONS 



A, ENGLISH LETTER SYMBOLS 

a Distance from pivot to the center of rotation (m) . 

A Area of one lining shoe (m^). 
b Width of friction material (m) • 

B. Biot modulus, 
c Specific heat (J/Kg-°C), 

C Thermal capacity (J/°C). 

d Distance from actuating force to the hinged pin (m) . 
dc Rata of deceleration (m/sec^ ) • 

E Kinetic energy (J) . 
f Frictional force (N) 

F Actuating force (N) • 

Fq Fourier modulus, 
g Gravity constant {m/secl ) , 

h Convection heat transfer coefficient (W/m^-C), 

k Thermal conductivity (W/m-°C), 

Friction moment (N-m) , 

M Normal moment (N-m) , 
n 

M Normal force (N| , 

p Pressure between lining and drum at any point (N/m^), 
p Maximum pressure between lining and drum (N/m^), 

B. 

Q Heat generated (W) , 
r Inside drum radius (m) , 

R Wheel radius (m) , 

Thermal resistance (°C/W), 
t Time (sec,) 
tk Thickness (m) , 

T Temperature (°C) , 

Tq Torque (N-m) , 

V Velocity (rn/ssc,), 

Vo Volume (m^), 

W Vehicle weight (N) , 
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B. NOTAriON 



R. . The thermal resistance between node i and the 
ij 

adjoining node j. 

tP The temperature of node i at time step p. 

C. GR2EK LETTER SYMBOLS 

9 The angle between the hinged pin and an element 
area on the lining, 

0 The angle at which the pressure between the lining 
and drum is maximum. 

Ii Friction coefficient, 
u Cold friction coefficient. 

Hot friction coefficient, 
a Thermal diffusivity (m^/sec.). 

A Finite increment. 
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I. 



INTRODUCTION 



Brakes are aiechanical de7ices for retarding the motion 
of a vehicle or machine by means of friction. Because of the 
similarity of their functions, many clutches may also be 
included here, assuming centrifugal forces are accounted 
for. 

A simplified dynamic representation of a brake is shown 
in Fig. 1. Two masses with inertias, and I 2 * rotating 
at the respective angaiar velocities o)i and 0)2 (one of which 
may be zero) , are to be brought to the same speed by 
engaging the brake. 

The friction brake has three basic elements; two 
opposing friction surfaces and a mechanism for forcing the 
friction surfaces into contact. Whenever a friction brake is 
engaged to join two members having relative motion, there is 
a period of slip which may last several seconds. This slip 
is one of the chief merits of the friction brake; it absorbs 
saocks and prevents excessive torsional stresses on the 
power transmission system. On the other hand, slip is the 
limiting factor in friction clutch and brake performance; 
for heat is generated in proportion to slip, torque 
transmitted, and period of slip. 

The following parameters are of interest in analyzing 
the performance of these devices; 

1. The actuating force. 

2. The torque transmitted. 

3. The temperature rise. 

4. The slip time. 

This report deals with Internal- Expanding Rim Brakes. 
This formulation also applies to internal-expanding clutches 
if centrifugal forces are accounted for. 
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II 



INTERNAL-EXPANDING RI.1 CLUTCHES AND BRAKES 



A. GENERAL MACHANICAL PRINCIPALS 

A brake or clutch assembly/ uses a brake shoe to which 
is attached a frictioQ material, called lining. The lining 
is riveted or bonded to the brake shoe as shown in Fig. 2. 
rhe brake shoe is pivoted at a fixed point and the other end 
is subjected to a force which presses the shoe in contact 
with the drum. The force between the brake and the drum is 
radial as the drum rotates. If a point on the rotating drum 
surface first makes contact with the shoe at the end nearest 
the pivot, the shoe is termed a "trailing shoe". If it first 
makes contact at the other end the shoe is termed "leading 
shoe", the latter giving a higher braking torque than the 
former for a given braking force. 

The friction between the lining and the drum creates 
heat which is basically the conversion of energy of motion 
of the vehicle or machine to thermal energy at the friction 
surfaces, namely the lining and the drum. This heat is then 
dissipated and absorbed by the drum by conduction, 
convection and radiation into the atmosphere. 

B. FRICTION FUNDAMENTALS AND MATERIALS 

Friction mechanisms, such as brakes, are systems for 
converting mechanical energy into heat. Several basic 
factors affect friction and wear of materials used in brake 
systems. The main factors are temperature, pressure, speed, 
surface roughness, and type of material. Some organic or 
molded friction materials show no change in friction 
characteristics with pressure, while others such as 
sintered- metal materials decrease in friction coefficient as 
pressure is increased. For metallic friction materials there 
is also a decrease in coefficient of friction as speed 
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increases. Temperatare effects upon the coefficient of 
friction vary widely with the type of materials used. 

In a two-shoe internal expanding brake there is a 
tendency for the brake drum to deform under hard 
application. Drums become elliptical and the force to do 
this is quite high and contributes to friction force. 

A brake or clutch friction material should have the 
following characteristics to a degree which is dependent 
upon the severity of the service; 

1. A high and uniform coefficient of friction. 

2. The ability to withstand high temperatures, together 
with good heat conductivity. 

3. Properties which are not affected by environmental 
conditions such as moisture. 

4. Good resiliency. 

5. High resistance to wear, scoring and galling. 

C. BRAKE DRUMS 

One of the primary functions of a brake drum is that of 
absorbing and dissipating the heat developed during the 
application of the brake. A brake drum is a heat sink into 
which heat goes after it is created by the rubbing friction 
of the brake lining contact to drum. The brake shoe and 
lining permanently fixed on the axle, when actuated, 
contacts the drum under pressure to cause the friction to 
stop the vehicle. The energy of motion of a vehicle is 
converted to thermal energy by the brake assemblies. A brake 
drum must have the capacity to absorb and dissipate this 
heat energy within the limits of the brake heat, input. If 
this is not the case, the drum expands and the brakes fade 
or fail. The greater the mass of the drum, the more heat it 
can absorb and store until such time as the heat can be 
dissipated by convection and radiation [Ref. 1]. 

An ideal brake drum would have the following 

characteristics ; 
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1. High structural strength to resist bursting forces. 

2. Uniforni coefficient of friction. 

3. Hard surface to resist scoring. 

4. High heat conductivity to rapidly conduct heat away 
from braking surfaces. 

5. High emissivity factor to radiate heat from the drum 
surface to the atmosphere. 

6. High heat storage capacity to store heat from 
successive brake applications until it can be 
dissipated. 

7. Good machinability to permit boring of the drum. 

D, STATIG AND DYNAHIG ANALYSIS 

1 . Assumptions 

In developing the equations^ the following 
assumptions have been made; 

a. The pressure at any point on the shoe is proportional 
to the moment arm of this point from the pivot. 

b. The effect of centrifugal force may be neglected. 

c. The shoe is assumed to be rigid. 

d. The friction coefficient is a linear function of 
temperature and it does not vary with pressure^ wear 
and environment. 

2. Pressure Concept 

To analyze an internal shoe refer to Fig. 2, which 
shows a shoe pivoted at a fixed point with the actuating 
force acting at the other end of the shoe. The mechanical 
arrangement does not permit pressure to be applied at the 
pivot, therefore the pressure at this point is zero. If the 
shoe rotates through a small angle about A, the radial 
movement of any point on the arc of contact, is proportional 
to the moment arm of this point from the pivot. Assuming 
that the material of the brake lining and support obey 
Hooke's law, the pressure at this point will also be 
proportional to this moment arm. The distance is 
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f f p 



proportional to sin 0 . Therefora, the relations between 
pressare at any point and the maximum pressure, p , will be 
given by the following formula; 

P ^ Pa 
sin6 sine 

a 

From this formula it can be seen that the frictional 
material at the heel, contributes very little to the braking 
action, therefore it is better to begin the friction 
material at an angle 0x greater than, say 0.15 rai. It can 
be seen also that the pressure will be maximum when 0 =90° 
or if the toe angle dz is less than 90°, then the pressure 
will be maximum at the toe. For good performance it is 
recommended to concentrate as much frictional material as 
possible in the neighborhood of the point of maximum 
pressure [Ref. 2]. 

3. Actuating Force and Torque Calculation 

From Fig. 2, it can be seen that the differential 
normal force on an element area of the lining will be; 

dN = pdA (2) 

where dA is an area element of the lining and it's magnitude 
is; 



dA = rbd0 



(3) 



In Equation 3, r is the inside drum radius and b is the drum 
width. Substituting for p and dA gives; 
p brsinB 



dN = 



sin0„ 



de 



(4) 



At the same point the differential frictional force is; 

df = lidN (5) 

where \x is the coefficient of friction. 
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h 




rhe actuating force, F, can be calculated using the fact 
that the summation of the moments about the hinge pin is 
zero. The moment due to frictional forces is; 



M- =q/ ^(r-cos9)df 

I 9 1 



( 6 ) 



where a is the distance from the pivot to the center of 
rotation. Substituting the value of dF and integrating 
from to 02 gives; 

a 

— {(COS0 1-COS02 )+ 2 p(sin^S i-sin^02 ) } 
a 

where y is assumed to be constant along the lining. 
Similarly the moment due to normal forces is given by; 



(7) 



,i =„/'^^asin 0 dN 
n D 1 



( 8 ) 



Substituting the value of dN and integrating 
from 01 to 02 gives; 
p bra 



M 



n Sint 



-{O. 5 ( 02 - 9 i)-O. 25 (sin 02 -sin 0 i) 



(9) 



The actuating force must balance the moments, therefore; 






( 10 ) 



where d is the distance from the hinge to the point of 
application of F, The torque applied to the drum by the 
brake shoe is; 



To =^/®^rdf 



( 11 ) 



After substituting the value of df and integrating ; 
UP„br^ 

To = — : 5 ( COS0 1 -COS0 2 ) 



( 12 ) 
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4. 



Rate of Heat Senerated and Deceleration Calculation 



The differential rate of heat generated by an 
element area of the lining is equal to the velocity of the 
inside surface of the drum relative to the lining , rimes 
the differential frictional force acting on the element 
ar ea ; 



dQ = V^df (13) 

Assuming the brake is on a vehicle wheel with a radius of R, 
the inside surface velocity is equal to; 




where 7 is the velocity of the vehicle and is a function of 
ti me. 

If 7 = 7 (t) then 7 =7 (t) and the heat generated will be also 
r r 

a function of time. Substituting the values of 7 and df and 

r 

integrating from 9 i to 02 r we get the following formula for 
the heat generated at any time t, 

Pgby 2 

Q(t) = sine (^)(cosSi-cos92)V(t) (15) 

a 

The kinetic energy of a vehicle of weight W is given by; 

E = |(|)V" (16) 

Mote that if the brake is on a four wheel vehicle, there 
will be eight shoes. Assuming all are leading shoes, each 
will stop one-eight of the vehicle weight, so W /8 must be 
used in Equation (16) . The rate of change in the kinetic 
energy is; 



dt 



dt 



(17) 



From the energy conservation law the rate of change in the 
kinetic energy is equal to the heat generated; 
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Q(t) 



dt 



(18) 



Substituting the valu 2 of Q (t) and dE/dt, it is seen that 
the velocity V (t) cancels and so the deceleration is not a 
function of tiae. Therefore the deceleration, dc, is; 

a 

The velocity at any time is; 



V = V.-dct 
1 



( 20 ) 



where is the initial velocity. Substituting the velocity 
in Equation (16) , yields the rate of heat generated as a 
function of time, 

P„bu 2 

Q(t) = -f-^(|)(COS0l-COS92)(V -dot) (21) 

Sin u ri 1 



In this study the friction coefficient was taken as constant 
up to a temperature of 90°C and after 90°C, decreases 
linearly to zero at a specified temperature, 



y = 



u -bi, 
h 

AT 



(T-90) 



T < 90 



90 °C < T < T 

— max , 



( 22 ) 



T > T 



where is the cold coefficient of friction and is the 

hot coefficient of friction. 



E. saRFASE TEMPERATURE CALCULATION 

Since the function of a brake is to convert kinetic 
energy into heat, surface temperatures of brake linings and 
drums are most important. Therefore it is necessary to know 
the temperature of the mechanism during and after any stop. 
The temperatures were calculated by the finite difference 
method. 
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1 . ^Assumptions 

a. One dimensional heat flow-The heat flow is from the 
inner surface to the outer surface of the drum. 

b. Constant heat transfer coefficient. 

c. No heat dissipated by radiation. 

d. The heat is generated on the inner surface. 

2. lemperature Analysis 
a. Theory 

The differential equation to be solved in order 
to find the temperature in the drum, based on the 
assumptions, is; 



3^T , Q ^ ,1.3T 
3x^' k 



(23) 



with the following boundary conditions: 
at x=0 heat is generated, 

at x=tlc heat is transfered to the atmosphere by 
convection. 

In the equation above k is the thermal conductivity, a is 
the thermal diffusivity, t is time and tk is the drum 
thickness. This ajuation can be solved by the finite 
difference method ]Sef. 3]. The finite difference model 
used here is shown in Fig. 3. The rate of change with time 
of the internal energy of a node i is approximated by; 



mPt ^_'T«P 

AE i i 

_ = pcAVo— 



(24) 



where p is the density, c is the specific heat and is 

the drum volume. 

Now define the thermal capacity as 



C. = p.o.AVo. (25) 

1 ^111 

The forward difference equation for all nodes and boundary 
conditions is; 
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« m 

m 










(26) 



qM 



rpP_^pP 



th, ij 



C.- 

1 



rpP+1 

i 



-Tr 

1 



At 



where 3^, . . is the thermal resisrance 

D+l 

Solving the above equation for gives; 



v^, . . C. 

th,ij 1 



1 th,ij 



)T^ 



(27) 



The thermal resistance can be calculated from the geometry 
and boundary conditions "Ref. 3]. To ensure 

stability At must be equal or less than the following 
nodal relation; 



At < ( 



C. 



(28) 



th, ij 

With the assumptions made, the drum can be viewed as an 
infinite plate, with heat generated at the surface of the 
first noie, as shown in Pig. 3. It is assumed that in every 
drum, there are two shoes and that both are leading shoes. 



Therefore, two times must be taken. 



h. (1- 

th,ij 1 



i ^th, ij 



■)T? 

1 



(29) 



b. Formulation 

In the computer program 5 nodes were taken. In 
order to oheck accuracy, the program was run with 7 and 10 
nodes. In each case the result was the same within 5 '^C. 
The heat is generated in the inner drum surface. Therefore Q 
appears in the formula of temperature in the first node and 
for all the other nodes Q is equal zero. With the 
assumptions mentioned above, the heat transfer through the 
drum is solved as a heat transfer problem through an 
infinite plate, with heat generation at the inner surface 
and with a heat convection boundary on the outer surface as 
shown in Fig. 3. Equation (29) can be simplified using two 
dimensionless parameters, Biot and Fourier moduli!. 
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B. 

1 



hAx 

k 



( 30 ) 



aAt 

(Ax)- 



(31) 



rhe final equations for calculating the temperatures at the 
nodes now become; 

For the first node; 

2qPa1 



rP+1 - 



+ (1-2Fo)tP + 2FoT| 



For the interior nodes; 

" 2)tP} 

For the last node; 

Tf 1 = 2 Fo{tP_^+B.T„+ (^ -B.-DtP} 



(32) 



(33) 



(34) 



F. BRAKE DUTY CYCLE 

In addition to tae parameters mentioned above the design 
of a brake depends on the initial speed, final speed, number 
of stops, and the rest time between each stop. In this 

analysis a general duty cycle was considered so that the 
initial speed, final speed and the acceleration period 
between stops can be different for each part of the design. 

In the design examples presented here, a vehicle was 



stopped 


four consecutive 


times with 


the follow: 




Initial 


Final 


Rest 




Speed 


Speed 






m/sec. 


m/sec. 


sec. 


1 


25.0 


0.0 


20.0 


2 


25.0 


0.0 


20.0 


3 


25. 0 


0.0 


20.0 


4 


25.0 


0.0 


- 
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OPTIi'ilZATION 



A. INTRODUCTION 

Enginsaring ana;i/sis using tha digital computer has 
become commonplace. It is less common to use the computer 
to ma!<a the actual design decisions, such as sizing of 
structural members or placement of mechanical linkages. This 
may be largely attributed to the fact that fully automated 
design requires techniques that are unfamiliar to much of 
the engineering community. 

In many engineering problems, it is necessary to 
determine the minimum or maximum of a function of several 
variables, limited by various linear and nonlinear 
inequality constraints. It is seldom possible, in practical 
applications, to solve these problems directly, and 
iterative methods are used to obtain the numerical solution. 
Machine calculation of this solution is, of course, 
desirable. The CONMIN program is available to solve a wide 
variety of such problems [Ref. 4]. 

CONMIN is a FORTRAN program, in subroutine form, for the 
minimization of a multi- variable function subject to a set 
of inequality constraints. The basic optimization algorithm 
is the Method of Feasible Directions [Ref. 5]. The user must 
provide a main calling program and an external routine to 
evaluate the objective and constraint functions and to 
provide gradient information. If analytic gradients of the 
objective or constraint functions are not available, this 
information is calculated by finite difference. While the 
program is intended primarily for efficient solution of 
constrained problems, unconstrained function minimization 
problems may also be solved, and the Conjugate Direction 
Method of Fletcher and Reeves is used for this purpose 
[Ref. 6]. 
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B. DEFINITION OF TERi^S 



Most iisciplines have a unique set of nomenclature used 
to describe the concepts within that discipline. Some of the 
commonly used terms in numerical optimization are summarized 
here. 

Objective- The value of the function which is to be 
minimized or maximized during the optimization process, 
synonyms are cost, merit and payoff. The common mathematical 
designation is F (X) . In the present study the objective was 
to minimize the material in the brake drum. 

Design variables- The parameters to be changed 
daring the optimization process in order to minimize or 
maximize the value of the objective function. Synonym; 
decision variables. The common mathematical designation is 
the vector X. Design variables considered in this study 
include, drum thickness, widthr the angle between the hinged 
pin and the end of the lining, and the distance from the 
pivot to the center of rotation. 

Inequality constraints- Dne-sided conditions which 
must be mathematically satisfied for the design to be 
acceptable. The common mathematical term is G(X)<0 or 
G(X)>0. If the ineqality condition is satisfied on G (X) , the 
design is acceptable, (feasible) . If it is not satisfied, 
the design is not acceptable (infeasible) . Constraints 
considered here include, vehicle stopping time, maximum drum 
temperature, and actuating force. 

Side constraints- Upper and lower bounds on the 
individual design variables X . The common mathematical 
representation is X^< X^< X^. 

Design space- The n-dimensional mathematical space 
spanned by the vector of design variables X. 

Active constraint- Constraint G.(X) is called active 
if its value is zero (or near zero for computational 
purposes) . 
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Inactive constraint- 

G (X)<0. 
j 

7iolated constraint- 

G (X)>0, 

j 



Constraint G . (X) 
J 

Constraint G . (X) 
J 



is inactive if 
is violated if 



C. THE OPTIMIZATION PROCESS 

The general design optimization problem can be stated 
mathematically as follows: Find the set of variables X^, 

1=1,2. ...n^ which will 

Minimize F(X ) (35) 



Subject to: 



G.(X) < 0 

< x“ 



j = l,2 m 

i=l ,2 n 



(36) 

(37) 



Vector X contains the set of independent design 
variables X i=l,2....n. X may represent, for example 
width, thickness, and angles in the brake optimization. The 
objective function used here is the drum volume. 

Equation (36) defines the ineqality constraints imposed 
on the design. For example, if the temperature on the inner 
drum surface must not exceed a specified value T, the 
associetei design constraint becomes, in normalized form 

li - 1 < 0 (38) 

fi 

The lower and upper bounds on the design variables, 
given by Eg. (37) , limit the region over which the functions 
F(X), and G (X) are defined. These constraints are often 
referred to as side constraints because they form the sides 
or bounds of the n-dimensional spaced spanned by the design 
variables X. 

If all the inequalities of Bqns. (36) and (37) are 
satisfied, the design is said to be feasible; if any of 
these conditions are not satisfied, the design is not 
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feasible. If F (X) is a miniaani and the design is feasible, 
it is alsD optimum, or at least, a relative optimum. Note 
that because the objective and constraints may be nonlinear, 
there may be multiple minima in the design space that cannot 
be identified using current methods. While this is a matter 
for concern, since it is desired to find the true optimum, 
it must be remembered that the same mathematical conditions 
exist if the design process is not automated. However, using 
optimization techniques, it is a simple matter to restart 
the optimization from several initial points in the design 
space and thereby improve the probability of obtaining the 
true optinum design, a process that would be quite time- 
consuming in manual design. 

Equations (35) -(37) define the nonlinear constrained 
optimization problem. If Eqs. (36) and (37) are not imposed 
on the design, the optimization problem is defined by 
Eq.(35) alone and is therefore an unconstrained minimization 
problem. 

Most nonlinear optimization algorithms update the vector 
of design variables by the iterative relationship; 

+ aS^ (39) 

where g is the iteration number, vector S is the direction 
of search in the design space, and the scalar a is referred 
to as a move parameter which , together with S, determines how 
much the vector X is changed during the g-th iteration. An 
initial design defined by X must be supplied. The 
optimization process then proceeds in two steps. First, 
the direction S, which improves the design, is found, and 
second, the scalar a , is determined which improves the 
design as much as possible when moving in this direction. 
The process is repeated until there is no further design 
improvement, indicating that this is the optimum attainable 
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1 



dssign. For further details see Ref. 7. 

D- COPES AND SUBROUTINE ANALIZ 

In order to siaplify the use of CONMIN and to further 
aid in the design optimization process a Control Program For 
Engineering Synthesis, COPES, was developed by Vanderplaats 
[Ref, 7]. COPES is the main program (recall that CONMIN is 
written in subroutine form) . The user must supply an 
analysis subroutine with the name ANALIZ, which will 
calculate the various parameters. This subroutine has three 
segments; INPUT, EXECUTION, OUTPUT. 

All parameters which may be design variables, objective 
functions or constraints are contained in a single labeled 
common block called SL03CN. 

Copes Terminology 

The COPES program currently provides six specific 
capabilities; 

1. Simple analysis, just as if COPES was not used. 

2. Optimization-Minimization or maximization of one 
calculated function with limirs imposed on other 
functions. 

3. Sensitivity analysis- The effect of changing one or 
more design variables on one or more calculated 
functions. 

4. Two-variable function space- Analysis for all specified 
combinations of two design variables. 

5. Optimum sensitivity- The same as sensitivity analysis 
except that, at each step, the design is optimized with 
respect to the independent design variables. 

6. Approximate optimization- Optimization using 
approximation techniques. Usually more efficient than 
standard optimization for up to 10 design variables or 
if multiple optimizations are to be performed [Ref. 7]. 
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IV 



DESCRIPTION OF TEE COMPUTER PROGRAM 



A. GENERAL PROGRAM ORGANIZATION 

A functional bloclc diagram o£ the program is presented 
in Fig, 4, A general description of the subroutines 
contained in the program is given here. Appendices A 
through D discuss the preparation of input data^ list the 
important computer program nomenclature^ and list the 

program. 

B. SaBROUTINES 

1 . Subroutine ANALIZ 

subroutine ANALIZ organizes the basic analysis used 
in the optimization. It controls the reading of the initial 
design description and calculation of the values of the 
objective function, constraints, and all other parameters 
necessary to solve the problem. COPES /CONMIN updates the 
design to minimize/maximize the objective function, 
iterating until no further improvement in the objective 

function is possible without violating one of the 

constraints. COPES/CDNMIN calls subroutine ANALIZ to obtain 
the function value luring the optimization. 

2. Subroutine INPUT 

This subroutine reads all input data associated with 
the brake analysis. Instructions for problem deck 

preparation are given in appendix B, 

3. Subroutine TEMPR 

This subroutine calculates the heat transfer 
constants such as the thermal capacity of each node and the 
resistance of each node, determines the time increment in 
order to insure a stable solution, and calculates the rate 
of heat generation. In order to calculate the temperature of 
each node, it calls two subroutines. From subroutine BEAK it 
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obtains the deceleration needed to calculate the rate of 
heat generated and from subroutine TEMA it obtains the 
temperature rise of each node. Then it calculates the 
temperatures during the time that the brake is not in use. 
This subroutine is also capable of calculating the 
temperature rise of a drum when a constant rate of heat 
dissipation is given. 

4. Subroutine rSdA 

This subroutine calculates the temperature of each node. 
As mentioned before^ the heat is generated on the inner 
surface^ and on the outer side of the drum the heat is 
dissipated by convection. The formulas used were developed 
by the finite difference method, and are given in section 

rr-s-2. 

5. Subroutine BRAK 

This subroutine calculates the torque, actuating force, 
and the friction moment of one shoe. It also calculates the 
drum volume and the deceleration of the machine. The 
subroutine takes into consideration a constant friction 
coefficient until a temperature of 90°C is reached and a 
linear decrease in the friction coefficient for higher 
temperatures. More details are given in section II-D-4. 

6. Subroutine paiPUT 

This subroutine echos the input data and prints out 
tne thermal and mechanical information for the brake. An 
example of the output obtained from this subroutine is shown 
in Table 1 and Table 2. 
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V. TEST PROBLEM AMD RESULTS 



The computer program was tested with the data specified 
in Table 1. The objective function which was minimized was 
taa volume of the drum material. Design variables were the 
drum width, the angle between the hinged pin and the end of 
the lining, the ratio of the pivot to center of rotation 
distance to drum radius, and the drum thickness. The side 
constraints (limitsi on the design variables were; 



Design Lower 

Variable Bound 

1. (3) Width, b 0.0 

2. (5) Theta 2, 1.2 rad. 

3. (12)Ratio,Rd 0.1 

4. (1 8) Thick ness , tk 40 mm. 



Upper 
Bound 
80 mm . 

2. 5 rad, 
0. 9 

No bound 



The number in parentheses is the location of 
variable in the COMMON block in the computer program. 
Constraints were imposed on the actuating force F, 



the 



the 



maximum temperature, T 

m 

and stopping time, t. 



on the inner surface of the drum 



Constrained 

Variable 

1. (9) Force, F 

2. (25) Time, t 

3. (4) Temperature,! 



Lower 
Bound 

200.0 N-m 
No bound 
No bound 



Upper 

Bound 

2500.0 N-m 

7.00 Sec. 

230.0 °C 



The vehicle which weights 25700.0 Newtons is stopped 
four consecutive times from a velocity of 90.0 Km/hr to 
zero, with an acceleration period of 20.0 sec. between 
stops. The values of the design variables and the 
constraints before and after optimization are; 
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Before After 

Optimization Optim izat ion 



Objective 
Function 
Drum Volume 


0.754 E-03 


0.159 S-^ 


Design 

variables 

Width 


0.08 m 


0.08 m 


Theta 2 


2.10 Rad. 


1.92 Rad 


a/r 


0.75 


0.755 


Thickness 


0.010 m 


0.020 m 


Constraints 

Actuating 

Force 


2815.1 N 


2086. 3 N 


Stopping 

time (last stop) 


7.04 sec. 


7.00 sec 


Temperature 
after last stop 


348.7 


229.20C 



Note that the objective function increased as a result 
of optimization. This is because the initial design violated 
constraints on stopping time and maximum temperature. 

Further results are listed in Tables 1 and 2. In 
addition to optimization^ a sensitivity analysis of the 
design variables and a two-variable function space analysis 
for width and thickness were performed. The graphical 
results are given in Figs. 5 through 17. The results can be 
summerized as follows; 

a. The affect of changing the inside drum radius with all 
other design variables held constant; 

As shown in Figs. 5-7, for small inside drum radii the 
drum temperature is very high. The stopping time is 
long and the torque is low. Inside drum radii over 130 
mm give reasonable drum temperature and stopping time, 
for the example considered. 
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b. As seen in Figs. 8-10, the affect of changing the druci 
ifidth with all other design variables held constant is 
the same as described above. 

c. The effect of changing the drum thickness with all 
other parameters held constant is; For a drum thickness 
up to 6 mm, the stopping time and drum temperature are 
considerably high. Over 16 mm thickness, the stopping 
tima re nains almost constant. For a small thickness the 
torgue is very low due to the high temperatures. For 
thicknesses over 20 mm, the torque remains about 
constant. 

d. The effect of changing the angle between the hinged pin 
and the end of the lining is; For a small 02 angle the 
stopping time is very long because the torque is low. 
The stopping tima becomes reasonable when 02 >1.8 Rad. 
Obviously there is an increase in the drum temperature 
as 02 increases but the overall change in temperature 
is small. 

e. From the two variable function space. Fig. 17, it can 
be seen that the constant volume line and the constant 
temperature line are almost parallel, this leads to the 
conclusion that for the cycle taken, the drum is a heat 
sink, and the amount of heat dissipated by convection 
during this cycle is small. 
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VI, TSflPSRATURS i^AISS 



SIi^PLIFIED CALCULATION 



A simplified way of finding the temperature rise of the 
drum is by using the equation; 

Q = |cAT (40) 

and setting Q equal to the amount of heat generated using 
Equation (21) from section II-D-4. This equation is in 
common use in engineering texts (See, for example. Refs. 2 
and 8) . The temperature rise calculated this way, is the 
average temperature of the drum, and not the temperature on 
the interface, which can be much higher (depending on the 
rate of heat generated) . Extreme temperature gradients cause 
distortion and excessive surface wear. Therefore it isn't 
always acceptable to use the simplified formula. From 
experience, it has been found that the surface wear 
increases dramatically as interface temperatures approach 
1*30 to 533 V (205 to 260 °C) , [Rsf. 8]. 

A comparison of the temperatures calculated on the inner 
drum surface, outer drum surface, and the average drum 
temperature calculated, using equation (40) , is given in 
Fig. 18. The graph shows the temperature rise for a vehicle 
stopped from a velocity of 90 xm/hr. From this graph, it can 
be seen that the drum temperature, based on equation (40), 
after the vehicle stopped is about the average temperature 
of the inner and outer surface temperatures. 

The results show that the drum will reach an uniform 
temperature of about 53°C, in 15 sec, after the vehicle has 
stopped. 

Calculating the temperature with the simplified formula, 
can lead to errors in the time needed to stop the vehicle. 
Because the temperature calculated with the simplified 
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fDrmala is lower than the temperature at the friction 
interface, the calculated friction coefficient is higher 
than the actual friction coefficient. Therefore the 
calculated stopping time will be shorter than the real 
stopping time. Ail this is true, provided the friction 
material behaves as assumed in section I-D-4. 

Because high temperature is detrimental to both the 
stopping ability and the wear characteristics of the brake, 
it is important that the interface temperature be calculated 
with reasonable accuracy in design. Fig. 18 clearly shows 
the temperature differences resulting from the two 
approaches . 

This difference in results is compounded when the 
simplified equation is used for design. Table 3 presents the 
design results based on the simplified approach. This design 
represents an apparent material savings of 27%. However, 
when this optimum is analysed using the finite difference 
heat transfer solution, the maximum temperature is 268.8 °C 
and the last stopping time is 7.54 sec. This time violates 
the constraint by 7.7%. Perhaps more importantly, the 
temperature at the interface of about 269 °C would surely 
lead to premature failure. Therefore this design is clearly 
too unconservative to be acceptable. 
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VII, CONCLUSION 



In sammary, a numerical optimization program is an 
effective way of finding a solution to an engineering 
problem, provided reasonable care is used in formulating the 
problem. 



yill. REC0.^M5NDATI0NS FOR FUTURE INVSSTIGATIONS 

The study has shown the feasibility of using numerical 
optimization in the design of Internal- Expanding Rim Brakes 
with two leading shoes. Further studies on the same design 
may be pursued by eliminating some of the restrictions. For 
example; 

1. To add heat dissipation fay radiation. 

2. To investigate drum temperatures for a drum with fins. 

3. To take into consideration changes in the surface 
pressure as a function of friction coefficient. 

4. To repeat all the calculations for a drum in which 
there is one trailing shoe and one leading shoe. 

5. To add the effect of centrifugal forces for clutches. 
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FINAL OPTIMIZATION INFORMATION ANA03670 
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APP2NDIX A 



LIST OF PARAMETERS 

A complete listing and description of all variables used 
in the program^ is not practical. The variables listed in 
this appendix are common to several subroutines of the 
program and will assist the reader in a study of the 
program. The Global location is the location of the 
parameter in the common block called GLOBCM, This common 
block is the means by which information is transfered 
between the subroutines and the CDPES/CONMIN program. 



Global 


Fortran 


Math, 


Location 


Na me 


Symbol 


1 


RI 


r 


2 


T 




3 


WDTH 


b 


4 


PRSA 


P 


5 


TETAl 


01 


6 


TETA2 


0 2 



7 


FRMNT 


Mf 


8 


ANMRT 


Mn 


9 


ACFRC 


F 


0 


n 


d 



Definition 



Inside drum radius (m) 
Torque of one shoe (N-m) 
Drum width (ra) 

Pressure between lining 
and drum (N/m^) 

The angle between the 
hinged pin and the (Rad,) 
begining of the lining 
The angle between the 
hinged pin and the end 
of the lining (Rad.) 
Friction moment (N-m) 
Normal momemt (N-m) 
Actuating force (N) 
Distance from actuating 
force to the hinged pin(m) 
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Heat generated (J/sec,) 
Distance from pivot to 
center of rotation (m) 
Cold friction coefficient 
Hot friction coefficient 
Friction coefficient at 
any temperature 
Drum surface area (m^) 
Outside drum radius (m) 
Drum thickness (m) 

An incremental thickness 
(m) 

Wheel radius (m) 

Car's weight (N) 
Deceleration (m/sec,^ ) 
Total time (sec.) 
Ecentricity (m) 
write statement control 
rime (sec. ) 

Drum volume (m^) 
Temperature at time p+1 
(sec. ) 

Write statement control 
Write statement control 
Write statement control 
Number of elements 
Number of segments 
Constant 

Gravitational constant 
Thermal conductivity 
(J/m-C) 

Convection heat coef f . 
(W/m2-C) . 

Specific heat (J/Kg-C) 
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TEMP 
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51-57 


RES 
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TC 
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FUR 
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66-72 
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VT 
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NELO 




81 


NELT 




82 


TETAA 
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83 


ACOF 




84 


TINT 
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85 


ZMAN 
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86 


NS HO 





Density {Kg . /m^ ) 

Time increment (sec.) 
ilax. temp, difference ( °C) 
Temp, at time p (sec.) 
Heat resistance ( °C/J) 
Heat capacity (J/°C) 

Biot moduli 
Fourier moduli 
Control parameter 
Control parameter 
Number of elements+1 
Number of elements+2 
The angle at which the 
pressure between the 
lining and drum is 
maximum. (Rad.) 

Constant 

Initial temperature ( °C) 
Time increment (sec.) 
Number of shoes 
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APPENDIX B 



INSTRUCTIONS FOR PROBLEM DATA PREPARATION 

Although the procedure is straight forward, preparation 
of input data for the program requires attention. Errors are 
easy to make and difficult to locate. Input data is 
described here for the brake analysis. For instructions on 
data preparation for optimization see Ref. 7. Input data 
should, in general, follow the steps outlined below. The use 
of the standard FORTRAN Eighty Columm Coding Sheet is 
recommended. Integer constants must be right justified in 
the appropriate field. Thera are eight input cards, read by 
subroutine INPUT, to describe the initial design, material 
properties and constants. Card format is given in 
parenthesis followed by specific instructions where 
necessary. 

1. First Card (110) - Duty cycle information. 

Cols 1-10 : Total number of consecutives stops 

and accelerations (NSEG) 

2. Second Card (I10,3F10.0) - Duty cycle information. 

a. Cols 1-10 : Control number. 

1 means-deceleration, 

2 means-brake not in use, 

b. Cols 1 1-20 : Velocity at start of deceleration. 

c. Cols 21-30 : The velocity at the end of the 

deceleration. 

d. Cols 31-40 : The time the brake is not in use. 

3. Third Card (5110) - Thermal analysis information. 

a. Cols 1-10 : Number of nodes (NSL) . 

b. Cols 11-20 : An integer number that controls the 

amount of printout when detailed 
output is required during the 
vehicle deceleration. The amount of 
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4. 



5 . 



c. Cols 21-30 



d. Cols 31-40 



e. Cols 41-50 

Fourth Card (7 

a. Cols 1-10 

b. Cols 11-20 

c. Cols 21-30 

d. Cols 31-40 



e. Cols 41-50 

f. Cols 51-60 

g. Cols 61-70 
Fifth Card (7F 
a. Cols 1-10 



lines written, depends on the 
stopping time and time increment. 
(NWR) . 

: An integer number that controls 
the amount of printout when 
detailed output of the temperatures 
is required during the period that 
the vehicle is not in use. The 
amount of lines written depends on 
the period length that the brakes 
are not used (NWRA) . 

: An integer number that controls the 
amount of printout when detailed 
output of the temperatures are 
required during the period of 
constant heat generation. (NWHQ) . 

: Number of braking shoes in the 
machine. 

F10.0) - Brake dimensions. 

: Inside drum radius (RI) . 

: Drum width (WDTH) . 

: Drum thickness (THK) . 

: Ratio of distance from pivot to 
center of rotation and inside 
radius (RD) . 

; Drum density (RHO) . 

: Angie between hinged pin and the 
begining of the lining (TETA1) . 

; Angle between hinged pin and the 
end of the lining (TETA2) . 

10.0) - Thermal and friction 
information. 

: Heat conduction coefficient (K) . 
(real number) . 
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b. Cols 11-20 : Heat convection coefficient (HCNV) . 

c. Cols 21-30 : Specific heat of the drum (SPHT) . 

d. Cols 31-40 : .dax, temperature difference between 

cold friction coefficient and hot 
friction coefficient (DFTH) . 

e. Cols 41-50 : Cold friction coefficient (CMIV) . 

f. Cols 51-60 : Hot friction coefficient (HMIV) . 

g. Cols 61-70 : Initial temperature (TINI) . 

6. Sixth Card (2F10.0) - Machine information. 

a. Cols 1-10 : Vehicles weight (W) . 

b. Cols 11-20 : Wheel radius (RTIER) . 

7. Seventh Card (5F10.0) - Analysis constants. 

a. Cols 1-10 : Maximum pressure between lining 

and drum (PRSA) . 

b. Cols 11-20 : Constant 3.1415927 

c. cols 21-30 : Gravitational constant (G) . 

d. Cols 31-40 : Increment of time (ZMAN) . 

e. Cols 41-50 ; The angle of maximum pressure 

(TETAA) . 

8. Sight Card (110. 0) - Print control. 

Cols 1-10 : An integer number can be zero or 1. 

If zero (or a blank card) - only 
the final results are printed. 

If 1- the temperature at time 
increments are printed. 
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STANDARD DECK STRUCTURE ANA05760 
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APPENDIX E 




Fig. 1 Dynamic Representation of a Brake or Clutch 
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Fig. 2 Brake Assembly 
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Drum Thickness 




Fig. 3 Finite Difference Model 
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Fig. 4 Block Diagram of the Program 
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Fig. 5 Maximum Inside Drum Temp. Vs. Inside Drum Radius 
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Stopping Time Vs. Inside Drum Radius 
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Fig. 7 Torque Vs. Inside Drum Radius 
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Fig. 8 Maximum Inside Drum Temp. Vs. Drum Width 
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Fig. 9 Stopping Time Vs. Drum Width 



1000r 




o 


C9 


Q 


C9 


Q 


GO 


Q 


O 


CO 


(D 




CVJ 



(w-N) snoaoi 



74 



Fig. 10 Torque Vs. Drum Width 
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Fig. 11 Maximum Inside Drum Temp. Vs. Drum Thickness 



16 h 




C03S) 3WI1 DNIddOlS 



76 



Stopping Time Vs. Drum Thickness 
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Fig. 13 Torque Vs. Drum Thickness 
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Fig. 14 Maximum Inside Drum Temp. Vs. Theta 
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Stopping Time Vs. Theta 
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Fig. 16 Torque Vs. Theta 
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Fig. 17 Two Variable Function Space 
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Fig. 18 Drum Temperature Vs. Time-Comparison 
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